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TERRESTRIAL PLANET FORMATION IN DISKS WITH VARYING SURFACE DENSITY PROFILES 

Sean N. Raymond', Thomas Quinn', & Jonathan I. Lunine- 
Short title: Terrestrial Accretion vs. Density Profile 

ABSTRACT 

The "minimum-mass solar nebula" (MMSN) model estimates the surface density distribution of the protoplane- 
tary disk by assuming the planets to have formed in situ. However, significant radial migration of the giant planets 
likely occurred in the Solar system, implying a distortion in the values derived by the MMSN method. The true 
density profiles of protoplanetary disks is therefore uncertain. Here we present results of simulations of late-stage 
terrestrial accretion, each starting from a disk of planetary embryos. We assume a power-law surface density 
profile that varies with heliocentric distance r as r"", and vary a between 1/2 and 5/2 (a = 3/2 for the MMSN 
model). We find that for steeper profiles (higher values of a), the terrestrial planets (i) are more numerous, (ii) 
form more quickly, (iii) form closer to the star, (iv) are more massive, (v) have higher iron contents, and (vi) have 
lower water contents. However, the possibility of forming potentially habitable planets does not appear to vary 
strongly with a. 

Subject headings: astrobiology — planets and satellites: formation — methods: n-body simulations 



1. INTRODUCTION 

\ In calculating the range of possible terrestrial planet systems 
. that could form around Sun-like stars, it is important to consider 
protoplanetary disks with surface density profiles different that 
the standard, "minimum mass solar nebula" model (MMSN - 
e.g., Weidenschilling 1977). Different disk profiles form pro- 
' toplanets with different masses and spacings (Kokubo & Ida 
. 2002). The accretion of terrestrial planets from these protoplan- 
' ets is strongly affected by the distribution of mass in the system, 
" and the dynamical timescales associated with regions of en- 
hanced density. The terrestrial planets that form in disks with 
varying density profiles may vary in significant ways, e.g. in 
their compositions and formation timescales. In this paper, we 
use a standard prescription for disk mass densities, constrained 
by observations of such disks, to explore how disk properties 
affect the growth of terrestrial planet systems. 

The MMSN model calculates the density profile of the pro- 
\ toplanetary disk in the following way. The mass in each planet 
i is augmented to Solar composition by adding the appropri- 
' ate amount of light elements. Next, this mass is smeared out 
into contiguous annuli. Finally, a simple, power-law profile 
is fit to the corresponding data points. The resulting surface 
density profile scales with heliocentric distance r (in AU) as 
E(r) = El r-3/2. Values for Ei range from 1700^ cm (Hayashi 
1981) to 4200 gcm-^ (Weidenschilhng 1977). In this case, Ei 
includes both the gas and dust component of the disk. The solid 
component of the MMSN follows the same r~^/^ profile with 
El « 6gcin~^. 

There is no reason to expect that all disks should follow 
the standard, r~^^^ MMSN profile. There is evidence that at 
least the outer planets in our Solar System may have under- 
gone significant radial migration (e.g. Levison & Morbidelli 
2003), which would distort the inferred density profile. In ad- 
dition, certain disk collapse models predict shallower, r"' pro- 
files (Shakura & Sunyaev 1973). In contrast, Kuchner (2004) 
showed that the mean, minimum-mass extra solar nebula has 
a steeper, roughly profile. Of course, this steeper profile 



also assumes that the known extra-solar giant planets formed 
and remained in situ. In light of studies of planet migration 
(e.g. Lin et at. 1996), this is not thought to be the case. There 
is clearly reason to consider protoplanetary disks with surface 
density profiles different that the MMSN model. 

Wetherill (1996) formed planets in both and r~^/^ pro- 
files, and saw no differences between the two cases. Chambers 
& Cassen (2002) simulated the late-stage accretion of terrestrial 
planets in two different surface density profiles. One of these 
was a simple, r~' profile, while the other had a low density in 
the inner disk, with a peak at roughly 2 AU, based on a temper- 
ature dependent condensation model of Cassen (2001). They 
found that the power law profile generally provided a better fit 
to the Solar system than the peaked profile, in terms of the num- 
ber and masses of terrestrial planets formed. There remains a 
great deal of uncertainty in disk models in terms of redistribu- 
tion of solids throughout the accretion process, including gas 
drag, type I migration, and dynamical scattering. 

Raymond, Quinn & Lunine (2004; hereafter RQL04) slightly 
varied the value of Ei, but only between E-IO g cm~^ . They did, 
however, vary the surface density in the outer protoplanetary 
disk by changing the planetesimal mass used in the outer ter- 
restrial zone. We found that a higher surface density results 
in the formation of a smaller number of less massive planets. 
These planets also tend to have higher orbital eccentricities. 
Raymond, Quinn & Lunine (2005; hereafter RQL05) ran two 
simulations with Ei = 2gcin~^. The systems that formed re- 
sembled extended asteroid belts, containing several Mars-sized 
bodies but no large terrestrial planets. 

We adopt a general power-law distribution E(r) = Ei 
(note that a bears no relation to the viscosity parameter by the 
same name). We consider three different power-law surface 
density distributions, with a = 0.5, 1.5 and 2.5. These cases 
encompass the likely range of actual values, according to the- 
ory and observation. We determine the effects of a on the ac- 
cretion process as well as the physical properties of the planets 
which form, including their mass, composition, and potential 
habitability. 
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2. INITIAL CONDITIONS 

We consider three surface density distributions. Each con- 
tains 10 Earth masses (M®) of material between 0.5 and 5 AU. 
Because of their different shapes, each has a different value for 
El. For a = 0.5, Si = 5.1 gcm"^, for a = 1.5, Si = I3.5gcm-^, 
and for a = 2.5, Si = 21.3^cm~^. The mass contained in the 
inner disk is much higher for higher values of a. The disk with 
a = 2.5 has 6.5 Mq of material inside 1.5 AU, compared with 
3.6 M0 for a = 1.5 and only 1.4 for a = 0.5. 

We assume in our initial conditions that oligarchic growth 
has taken place, and the disk mass is dominated by planetary 
embryos, shown to be the case for different values of a by 
Kokubo & Ida (2002). Each planetary embryo has swept up 
the mass in its feeding zone, which has a width of A mutual 
Hill radii. With our generalized power-law density profile, the 
mass of a given planetary embryo scales with its semimajor axis 
a and the separation between embryos A as 

M,„,,,„CXflP-")3/2^3/2j.3/2^ 

For a density profile steeper than a = 2, the mass of planetary 
embryos decreases with heliocentric distance, and the mass is 
very concentrated close to the star Conversely, shallower den- 
sity profiles result in higher-mass planetary embryos at large 
distances, with more mass in the outer disk. Figure 1 shows 
the distribution of evenly spaced planetary embryos for each of 
our surface density distributions. For a = 0.5, embryos are very 
small in the inner disk and quite massive (up to M®) and well- 
spaced in the outer disk. Conversely, for a = 2.5, embryos are 
smaller in the outer disk, though better- separated, as the Hill 
radius has a stronger dependence on semimajor axis a than on 
mass M (Rh ck a 

In Fig. 1, A is fixed at 7 Hill radii, to demonstrate the differ- 
ences in embryo mass and spacing. However, in our simulations 
we space embryos randomly by 1-2 Hill radii from 0.5 to 5 AU. 
The number of embryos varies for each value of a: we include 
roughly 380 embryos for a = 0.5, 290 embryos for a= 1.5, and 
280 embryos for a = 2.5. For each value of a, we perform three 
simulations with different random number seeds. 

The starting water contents of embryos are as described in 
RQL04: embryos are dry inside 2 AU, contain 0.1% water be- 
tween 2-2.5 AU, and 5% water past 2.5 AU. Their initial iron 
contents are interpolated between the values for the Solar Sys- 
tem planets and chondritic meteorites associated with known 
classes of primitive asteroids (values from Lodders & Fegley, Jr. 
1999). Embryos are given small initial eccentricities (e < 0.02) 
and inclinations (i < 1°). 

In each case we include a Jupiter-mass giant planet at 5.5 AU 
on a circular orbit. One could argue that if the mass of a giant 
planet scales with the isolation mass for a given surface density, 
then the masses of the giant planets should vary in position and 
mass according to our values of a (Lissauer 1995). However, 
we do not know which profile to be the most common in nature, 
or whether there even exists a universal density profile for pro- 
toplanetary disks. We are constrained by observations of actual 
extra-solar planets which do form, whatever their initial disk 
profile. Due to the large amount of uncertainty in this issue, it 
seems reasonable to assign the same outer giant planet to each 
disk. 

Each simulation was integrated for 200 Myr using the hybrid 
algorithm of the Mercury integrator (Chambers 1999), with a 6 
day timestep. Collisions are treated as inelastic mergers which 
conserve linear momentum, mass, water and iron content. Each 



simulation took 1-3 months to run on a desktop PC, and con- 
served energy to at least one part in 10^. 

3. RESULTS 
3.1. Details of two Simulations 

Figures 2 and 3 show the time evolution of two simulations 
with a = 0.5 and 2.5, respectively. The overall progression of 
both simulations is similar, but the details are quite different. 
In both cases, the eccentricities of planetary embryos increase 
until their orbits cross and collisions may occur. This eccentric- 
ity pumping occurs due to secular and resonant perturbations 
from both other embryos and the giant planet at 5.5 AU (which, 
for simplicity, we now refer to as Jupiter). Many embryos are 
ejected after close encounters with Jupiter. Larger bodies grow 
by accreting smaller ones, on timescales which depend on both 
the dynamical timescale (orbital distance) and the local density 
of material. Many bodies which form inside 2 AU contain one 
or more water-rich embryos which began the simulation past 2- 
2.5 AU. In time, most small bodies are removed from the sim- 
ulation, either via collisions with the forming planets, the Sun, 
or Jupiter, or via ejection from the system. What remains are 
a few (typically 2-4) terrestrial planets with significant masses, 
generally inside 2 AU. The details of the planets which have 
formed in Figs. 2 and 3 are summarized in the captions. 

The main difference between the two simulations has to do 
with their different mass distributions, and can be easily seen 
in a comparison of the eccentricities of the inner disk in the 0. 1 
Myr panels of Figs. 2 and 3. In 2, where a = 0.5, the majority 
of the terrestrial mass is in the outer disk, past 2.5 AU. Eccen- 
tricity pumping occurs mainly in the outer disk, via both mu- 
tual encounters and interactions with Jupiter, and the terrestrial 
planets therefore form from the outside in. The planets which 
form contain a large amount of water-rich material, simply be- 
cause so much initial material was found in the outer disk. The 
final planets contain only a small fraction of the initial mass, 
because embryos in the outer disk are very strongly excited by 
Jupiter, and largely ejected from the system. 

In Fig. 3, where a = 2.5, the majority of the mass lies in- 
side 2.5 AU. The eccentricities of embryos in the inner disk are 
self-excited through mutual gravitational interactions. Accre- 
tion proceeds quickly in the inner regions due to both the high 
density of material and short dynamical timescales. Terrestrial 
planets therefore form from the inside out. Most of the mass 
in the system is well-separated from Jupiter, so the final planets 
contain the majority of the initial mass. However, since there 
was little initial water-rich mass past 2.5 AU, they are much 
drier than in Fig. 2. 

3.2. Trends with a 

Figure 4 shows the final configuration of all nine simulations 
we performed, with the Solar System shown for comparison. 
Our simulations with a= 1.5 do not reproduce the Solar Sys- 
tem planets as well as previous simulations with the same pro- 
file, such as RQL04. This is because the total mass in the disk is 
higher in these simulations. The 10 M0 in the terrestrial region 
is roughly twice as much as in the MMSN model, and ~50% 
higher than in RQL04 (higher still for the case (ii) simulations). 
As shown in RQL04, a higher surface density of material re- 
sults in the formation of a larger number of more massive plan- 
ets with larger eccentricities. We see the results of this effect 
in comparing the a = 1.5 simulations with the Solar System in 
Fig. 4. None of our simulations form ~Mars-sized planets in- 
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side 2 AU, also because of the large total mass in protoplanets. 
Although this affects certain comparisons between our results 
and the Solar System planets, it does not affect comparisons 
between simulations with similar initial conditions. 

There exist several strong correlations between the final plan- 
etary systems and the slope of the density profile, a, which are 
shown in Table 1 . Perhaps the most striking is simply the lo- 
cation of the planets formed. In cases with a = 1.5 and 2.5, 
the innermost planet typically resides around 0.5 AU, while the 
innermost planet in the a = 0.5 simulations is usually around 
0.75 AU. It is important to note that in this analysis we define a 
planet to be a body larger than 0.2 Mq whose orbit lies inside 
2 AU. 

The trends from Table 1 can be summarized as follows. As 
compared with a protoplanetary disk with a shallow surface 
density profile, a disk with a steeper surface density profile 
forms a larger number of more massive terrestrial planets in a 
shorter time. Figure 5 shows the mean time for planets to reach 
50%, 75%, and 90% of their final mass for all three values of 
a. The reason that a = 0.5 planets form more slowly is that 
the timescale for accretion scales with the dynamical time and 
the local density. The density in the outer disk is high enough 
that planets form more from the outside in. The long dynam- 
ical timescales translate to slower planet growth, and planets 
at larger heliocentric distances. In the case of a = 2.5 planets, 
the high inner density and fast dynamical timescales conspire 
to form planets on remarkably short timescales. There is a con- 
straint from measured Hf-W ratios that both the Moon and the 
Earth's core were formed by f « 30 Myr (Kleine et al. , 2002; 
Yin et al. , 2002), suggesting that the Earth was at least 50% of 
its final mass by that time. Fig. 5 shows that planets in each of 
our density profiles satisfy this constraint. 

A larger number of planets is formed in the simulations with 
a = 2.5 because the planets in those simulations occupy a larger 
range in semimajor axes, i.e. they extend to smaller orbital radii 
than for a = 0.5. The mean orbital radii of the innermost planet 
for a = 2.5 and a = 0.5 simulations are 0.45 AU and 0.74 AU, 
respectively. This corresponds to a difference in orbital period 
of over a factor of two, providing enough dynamical room for 
another planet. Indeed, the a = 2.5 simulations formed an av- 
erage of three planets per simulation, as compared with two 
planets per simulation for a = 0.5. 

The reason that more mass ends up in planets for steeper den- 
sity profiles is the dynamical presence of Jupiter. During accre- 
tion, a large fraction of the material close to Jupiter is ejected 
from the planetary system. The "ejection range" of Jupiter, the 
zone from which a large fraction of material is ejected, extends 
to a factor of 4-6 in orbital period (a factor of 2.5-3.3 in semi- 
major axis) from Jupiter. If we naively assume as in RQL05 
that all material within a factor F in orbital period of Jupiter is 
ejected from the system, then the total mass ejected is 



27rl]i 



(2) 



where rj is the orbital radius of Jupiter. The value of F depends 
on the mass and eccentricity of the giant planet, and may also 
scale with rj and a. If we only consider material ejected from 
interior to the giant planet, then 



Mejec = n " ( 1 

2-a 



^l_p2/3ia-2)^ _ 



(3) 



We know the values of Ei and the mean amount of mass 
ejected per simulation for each value of a, so we can calculate 



a value for F (we also take into account the fact that embryos 
exist only out to 5 AU, while r/ = 5.5 AU). We find that F ranges 
between roughly 4 (a = 0.5) and 6 (a = 2.5). The mean loca- 
tion of the outermost planet in each simulation, independent of 
a, was at about the 5:1 resonance with Jupiter (at 1.9 AU), con- 
sistent with our values of F. The range at which Jupiter ejects 
material from an a = 2.5 disk appears to be farther than for an a 
= 0.5 disk (Jupiter's effective "ejection range" is the 6: 1 vs 4: 1 
resonance for a = 2.5 and 0.5, respectively). The reason for this 
increased range for a = 2.5 is the lack of material in the outer 
disk with which excited embryos may interact, possibly damp 
their eccentricities, and remain in the disk. A particle excited by 
Jupiter in an a = 0.5 disk may have its eccentricity subsequently 
decreased via embryo-embryo interactions, whereas the dearth 
of material in the outer disk prevents this from happening for a 
= 2.5. 

We perform a MMSN-like test on the planets which formed. 
For each system, we select all planets of 0.2 or more and 
spread their mass into concentric, nested annuli, and calculate 
an effective surface density. We choose the division of annuli 
between two planets to be the geometric mean of the semimajor 
axes of the two planets, following Kuchner (2004). For a = 2.5, 
we find a best-fit slope to the surface density, a fit, of almost ex- 
actly 2.5. However, we obtain a/u values between 2.1 and 2.3 
for the planets formed in a = 0.5 or 1 .5 simulations. The rea- 
son for this is because, as discussed above, Jupiter's ejection of 
material from the asteroid region has a much stronger effect on 
shallower density profiles. This decreases the surface density 
measurements for outer radii, skewing the fit to a steeper pro- 
file. If we only include planets inside 2 AU, outside of Jupiter's 
ejection range, the a/u for the a = 1.5 planets changes to 1.53. 
However, a/n for the a = 0.5 planets remains above 2, because 
of Jupiter's pervasive effects in low-a simulations. Choosing 
an even smaller outer radius would hkely yield a value closer to 
0.5, but the number of planets is too small to perform the test. 
This analysis implies that, in the absence of external pertur- 
bations (e.g. from the giant planets), terrestrial accretion forms 
planets that conform to their initial density profile. In this sense, 
the planets "remember" the initial state of the disk. 

This begs the question of why the Solar system planets fol- 
low an r"'-^ profile. Two explanations are possible: 1) the sur- 
face density profile of the protoplanetary disk followed an r"'-^ 
profile, or 2) the initial profile was different, but perturbations 
caused the final planets to follow an altered profile. If perturba- 
tions did affect the final planets, then the starting mass profile 
was likely flatter, i.e., the starting value of a was smaller than 
1.5. However, as demonstrated above, the ejection range of 
Jupiter extends inward to orbital radii of ~ 2 AU. Additional 
secular perturbations arise from the presence of Saturn (e.g. the 
secular resonance at 2.1 AU). But the inner terrestrial re- 
gion remains relatively sheltered from these perturbations. In 
addition, the giant planets' eccentricities may have been signif- 
icantly lower during the epoch of terrestrial accretion (Tsiganis 
et al. 2005), reducing the effect of secular resonances. We con- 
clude that the inner solar nebula probably did follow an r"'^ 
profile. 

The probability of a habitable planet forming between 0.8- 
1.5 AU with water mass fraction > 10~^ does not depend 
strongly on the value of a, although planets that form in disks 
with different a values will certainly have vastly different char- 
acteristics. Planets that form with a = 2.5 have much higher 
iron contents than those in a = 0.5 disks, simply because the 
amount of iron-rich material close to the star is much higher. 
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Table 1 
Accretion Trends with a 



Ircncl 


— n s 


(JL — V .J 


Ct — L.J 


Mean innermost planet' 


0.74 AU 


0.47 AU 


0.45 AU 


Mean planet mass 


1.4 Me 


1.7 Me 


2.2 Me 


Mean formation time^ 


55 Myr 


42 Myr 


18 Myr 


Mass (planets) / Mass (disk)-' 


0.28 


0.44 


0.65 


Total mass in planets 


2.8 Me 


4.4 Me 


6.5 Me 


Mean number of planets 


2.0 


2.67 


3.0 


Mean number of planets < 1 AU 


1.0 


1.3 


2.0 


Avg. mass of planets < 1 AU 


1.65 Me 


2.06 Me 


2.24 Me 


Mean largest planet 


1.75 Me 


2.8 Me 


3.0 Me 


Mean orbital eccentricity 


0.07 


0.12 


0.06 


Mean water mass fraction 


2.6 X 10"^ 


8.1 X 10-3 


1.8 X 10-3 


Mean iron mass fraction 


0.19 


0.27 


0.29 


Avg. num. of habitable planets'' 


1.0 


2/3 


2/3 



'We define a planets to be more massive than 0.2 Earth masses and to have semimajor axis a < 2 AU. 
•^Mean time to reach 90% of a planet's final mass. See Fig. 5. 
^Fraction of initial disk mass contained in final planets 

''Habitable planets are defined to have semimajor axes between 0.8 and 1.5 AU, and water mass fractions larger than 10 



Similarly, planets that form in a = 0.5 disks have much higher 
water contents than those in a = 1.5 or 2.5 disks, because the 
bulk of the solid mass in a = 0.5 disks is found in the water- 
rich outer regions. The iron and water contents will have large 
consequences for the nature of the surface and interior of these 
planets, as well as their potential habitabihty. 

4. DISCUSSION 

The true surface density profiles of the solid component of 
protoplanetary disks are unclear. Their formation involves a 
great deal of physics, including the condensation sequence as 
a function of temperature, and radial migration via gas drag 
and type I migration. It is hkely that the standard, minimum 
mass solar nebula model is oversimphfied, as there is strong 
evidence from the dynamical structure of the Kuiper belt that 
the giant planets experienced significant radial migration (Lev- 
ison & MorbidelU 2003). 

We have demonstrated that the terrestrial planets which form 
in protoplanetary disks with three different power-law surface 
density profiles have very different characteristics. The total 
number and mass of terrestrial planets is strongly affected by 
the surface density slope, a. In the presence of an outer giant 
planet, a disk with a steeper profile (higher a) yields a larger 
number of more massive planets than a shallower profile. Plan- 
ets in steeper-profiled disks form more quickly, have higher 
iron content and lower water content than for shallower pro- 
files, with imphcations for the habitability of planets in each 
case. 

It is unclear whether our placement of a Jupiter-mass giant 
planet at 5.5 AU is self-consistent for each chosen value of a. 
Kokubo & Ida (2002) and Lissauer (1995) argued that the char- 
acter of the giant planets may depend on the disk's density pro- 



file. We chose identical initial conditions to explore the effects 
of a in a "typical" planetary system, although the nature of such 
a system may be inextricably linked to a. Some of our results 
that were strongly influenced by the presence of a giant planet, 
such as the strong depletion of the a = 0.5 disk, may only hold 
under our given assumption of a giant planet at 5.5 AU. How- 
ever, many of our results are independent of this assumption, 
such as the formation timescales and compositions of the plan- 
ets we form. 

Using previous results, we can extrapolate the effects of other 
parameters on our simulations, such as giant planet orbital ra- 
dius and eccentricity, and the location of the snow line. It has 
been shown that a giant planet on an eccentric orbit forms ter- 
restrial planets which are drier, less massive, and have larger 
eccentricities than for a giant planet on a circular orbit (Cham- 
bers & Cassen, 2002; RQL04). In the presence of an eccen- 
tric giant planet, more of the material from the asteroid region 
is ejected, and planets form at a larger distance from the giant. 
Therefore, an eccentric giant planet would cause even more dis- 
parity between disks with shallow vs. steep density profiles. A 
disk with a = 0.5 would have an even higher fraction of its total 
mass ejected, causing the terrestrial planets to be smaller. An 
a = 2.5 disk, however, would not be strongly affected, since its 
mass is so centrally concentrated. The primary consequence of 
the giant planets' orbital radius is the amount of sohd material 
destroyed through ejection from the system. This will vary with 
heliocentric distance r and a as shown in equation 2. The loca- 
tion of the snow fine determines how much water-rich material 
exists in disks with different values of a. 

Kuchner (2004) found that applying the MMSN technique to 
the known, extra-solar planets yielded a very steep density pro- 
file with a ~ 2. This steep profile corresponds to an in situ for- 
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mation scenario for close-in giant planets (e.g. "hot jupiters"). 
If such a profile is common, then perhaps hot jupiters exist in 
these systems. We expect that terrestrial planets may form in 
the presence of a hot jupiter, even if it has migrated through the 
terrestrial region (RQL05). Indeed, terrestrial planets appear 
to be able to form in some of the known extra solar planetary 
systems, such as 55 Cancri (Raymond & Barnes 2005). 

Real protoplanetary disks are not perfectly smooth power 
laws. The physics and chemistry determining their true shape 
is very complex, and has not been fully modeled. Some pre- 
liminary models show a peak in surface density around 1-2 AU 
(Cassen 2001), which may explain the large masses of Venus 
and Earth (Chambers & Cassen 2002). Such a profile can be 



roughly approximated as having a small a value in the inner 
disk, and a larger a value past ^ 2 AU. We have shown that ac- 
cretion proceeds from the inner disk outwards for large a and 
from the outer disk inwards for smaller a values. In this com- 
posite disk, accretion would likely start in the highest density 
region around 1-2 AU, and proceed both inward and outward. 
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Fig. 1 . — Distribution of planetary embryos for our three surface density profiles. In this plot, embryos are spaced evenly by 7 mutual Hill radii (i.e., A is fixed at 
7). In our actual simulations, A varies randomly between 1 and 2, increasing the number and decreasing the individual masses of embryos, as compared with this 
figure. 
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Fig. 2. — Six .snapshots in time of tlie evolution of a simulation with a = 0.5, which formed two planets inside 2 AU. The size of each body represents its relative 
physical size, but is not to scale. The color represents its water content, and the dark circle in the center represents the size of its iron core. The inner one has 
semimajor axis a = 0.79 AU, eccentricity e = 0.08, mass M = \ A M^, and a water mass fraction W.M.F. of 1.6x 10"-. The second planet has a = 1.38 AU, e = 0.04, 
M=1.7 MfB, and W.M.F. = 2.7 xlO^^. 
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Fig. 3. — Six snapshots in time of the evolution of a simulation with a = 2.5, which formed three planets inside 2 AU, formatted as in Fig. 2. The inner planet has 
a = 0.52 AU, e = 0.03, M = 3.4 M®, and W.M.F. = 2.3 X Ur'. The second has a = 0.94 AU, e = 0.03, M = 0.8 M®, and W.M.f. = 3.1 X lO^^. The third planet has 
a= 1.40 AU, e = 0.05, M= 1.8 Mg,, and W.M.F. = 3.7 xlQ-^. 
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